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ABSTRACT 

We discuss MAXIPOL, a bolometric balloon-borne experiment designed to measure the E- 
mode polarization of the cosmic microwave background radiation (CMB). MAXIPOL is the 
first bolometric CMB experiment to observe the sky using rapid polarization modulation. To 
build MAXIPOL. the CMB temperature anisotropy experiment MAXIMA was retrofitted with a 
rotating half-wave plate and a stationary analyzer. We describe the instrument, the observations, 
the calibration and the reduction of data collected with twelve polarimeters operating at 140 GHz 
and with a FWHM beam size of 10 arcmin. We present maps of the Q and U Stokes parameters 
of an 8 dcg^ region of the sky near the star /3-UMi. The power spectra computed from these maps 
give weak evidence for an EE signal. The maximum- hkelihood amphtude of £{£ + l)Cf^/2 IT IS 
55II5 /xK^ (68%), and the likelihood function is asymmetric and skewed positive such that with a 
uniform prior the probability that the amplitude is positive is 96%. This result is consistent with 
the expected concordance ACDM amplitude of 14 /iK^. The maximum likelihood amplitudes 
for e{i+l)Cf^/2Tr and £{£ + l)Cf^ /2Tr are -'SltH and 18^3^ /iK^ (68%), respectively, which 
are consistent with zero. All of the results are for one bin in the range 151 < £ < 693. Tests 
revealed no residual systematic errors in the time or map domain. A comprehensive discussion 
of the analysis of the data is presented in a companion paper. 

Subject headings: CMB polarization, polarimetry, half-wave plate 
1. Introduction 



^School of Physics and Astronomy, University of Min- 
nesota, Minneapolis, MN, 55455, USA 

^Department of Physics, University of CaUfornia, Berke- 
ley, CA, 94720, USA 

^School of Physics and Astronomy, Cardiff University, 
Cardiff, UK, CF24 3YB 

^Jet Propulsion Laboratory, Pasadena, CA, 91109, USA 

^Astrophysics, University of Oxford, Oxford, UK, OXl 
3RH 

® Physics Division, Lawrence Berkeley National Lab, 
Berkeley, CA, 94720, USA 

'^Computational Research Division, Lawrence Berkeley 
National Lab, Berkeley CA, 94720, USA 

^Astrophysics Group, Blackett Lab, Imperial College, 
London, UK, SW7 2AZ 

^Space Sciences Laboratory, University of California, 
Berkeley, CA, 94720, USA 

'^''Department of Particle Physics, Weizmann Institute of 



Measurements of the polarization of the cos- 
mic microwave background radiation (CMB) can 
confirm fundamental predictions made by our cos- 
mological model and probe the period after the 
Big Bang when inflation is believed to have oc- 
curred. Several experiments recently reported 
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detect ion of either the TE or EE power spec- 
trum (iReadhead et al.ll2004l: iBarkats et al.ll2005t 



Montrov et al.ll2005HLeitch et al.ll2005HPage eral 



20061 ) and there is on-going effort to build even 
more sensitive receivers that can probe the small 
signal amplitudes that are predicted for the BB 
power spectrum. 

MAXIPOL is a bolomctric. balloon-borne ex- 
periment designed to measure the polarization 
of the CMB. It uses a rotating half- wave plate 
(HWP) and a wire grid analyzer to modulate 
signals reaching the bolometers from polarized 
sources. MAXIPOL is unique in that it is the 
only bolomctric CMB experiment to date to de- 
ploy rapid polarization modulation. With our 
technique, the polarized component of the inci- 
dent radiation is modulated at a frequency equal 
to four times the rotation frequency of the HWP. 
This property is advantageous for measuring the 
small CMB polarization because polarized sky sig- 
nals are moved to a narrow band in the frequency 
domain which can be far from detector 1// noise 
or other spurious instrumental signals. In partic- 



ular, there is no confusion with spurious signals 
that appear at the modulation frequency itself. 
Because only polarized signals are modulated, a 
HWP polarimctcr separates the temperature and 
polarization signals in the frequency domain, mak- 
ing analysis of the signals independent. Of partic- 
ular importance is the property that each detec- 
tor is an independent polarimeter measuring the 
Stokes parameters /, Q, and U over a relatively 
short time. Several systematic errors that compli- 
cate detector differencing techniques are avoided. 

Differencing polarimeters are defined here as in- 
struments that measure Stokes parameters by dif- 
ferencing the signal from two detectors, each of 
which is sensitive to one of the orthogonal lin- 
ear polarizations. Differencing polarimeters can 
be constructed using individual polarizers placed 
at the entrance aperture of each photometer, two 
bolometers with orthogonal absorbing grids, or or- 
thomode transducers which split the incoming po- 
larization into two components (see for example 
Masi et all l2006t iBenoit et all 12004) . Differenc- 



ing polarimeters only modulate polarized sky sig- 
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nals through sky rotation and telescope scanning. 
They are prone to spurious polarization signals 
through errors in the absolute calibration of detec- 
tor pairs, time dependent responsivity variations, 
noise properties that are not common to both de- 
tectors, or different antenna patterns for the two 
detectors. Polarization modulators like the one 
used in MAXIPOL can mitigate these problems. 

Polarization modulation technologies differ- 
ent from MAXIPOL's include rotating polariz- 
ers, photoelastic modulators and Faraday rota- 
tion modulators. Rotating polarizers reflect the 
unwanted polarization component and can thus 
produce spurious signals through multiple reflec- 
tions inside the receiver. Suitable photoelastic 
modulators have not been developed in the fre- 
quency bands of interest. Broadband Faraday 
rotation modulators for the mil limeter-wave band 
are just now being developed ( Yoon et al. 20061 ) 
and little is known about systematic errors associ- 
ated with their operation. Modulation of sky sig- 
nals with a rotating HWP and stationary analyzer 
is a proven techni que in in frared and millimeter- 
wave astrophysics (|Tinbergeniil996,) and therefore 
we chose to implement the technique for CMB 
polarimctry Since MAXIPOL is the first CMB 
experiment to produce results using this strategy, 
the experience gained from hardware implementa- 
tion, data analysis and characterization of HWP- 
specific systematic errors will inform the design of 
next-generation experiments that aim to charac- 
terize the anticipated B-mode signals. 

To build MAXIPOL, the receiver from the 
CMB temperature anisotropy experiment MAX- 



IMA (iHananv et al.l20 00':'Balb i et al.l2000l:lLee et al 



2001 



Stompor et al.ll2001. : .Abroe et al.ll2004[ ) was 



converted into a polarimeter by retrofitting it with 
a HWP and a fixed polarization analyzer. The rest 
of the MAXIMA instrument including the detec- 
tor system, cryogenics, optics, and electronics was 
essentially unchanged. MAXIPOL observed the 
sky during two flights that were launched from 
the NASA Columbia Scientific Ballooning Facility 
(CSBF) in Ft. Sumner, New Mexico. The first 
flight, launched in 2002 September, yielded less 
than an hour of useful data because of a telemetry 
failure. The primary CMB data set, which will be 
discussed in this paper, was collected during the 
second flight in 2003 May. 

We describe the instrument and the observa- 



tions in Sections [2] and [3l The HWP-specific in- 
formation appears in Sections [H [5l [HI and [51 In 
particular. Section [5] describes the processing of 
our time-domain data. This information will be 
useful for future CMB experiments that will use 
HWP polarimeters. Estimated power spectra and 
Q and U maps are presented in a summary of our 
analysis in Section[7l A co mprehensive des cription 



of the analysis is given in Wu et al. ( 2006h . 
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Fig. 2. — A cross-sectional view of the receiver near 
the focal plane showing the polarimetry components. The 
HWP was mounted at an aperture stop of the optical sys- 
tem and a roof-shaped polarization analyzer was mounted 
in front of the entrance apertures of the horns. A sketch 
showing the analyzer and the focal plane projected on a 
perpendicular plane is given in Figure [3] 



2. Instrument 

In this Section, we give an overview of the ex- 
periment emphasizing the polarimeter. Additional 



Dublications (Lee et al. 1998; Hananv et al.l 


2000l 


Rabii 20021 Winant 


20031: Rabii et al. 2006 


) and 


MAXIPOL papers ( 


Johnson et al.l 20031 Johnson 


2004; Johnson et al. 


20061 Collins l2006h. 
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MAXIPOL employed a three-mirror telescope 
with a 1.3 m ofF-axis parabolic primary mirror. 
The elliptical secondary and tertiary reimaging 
mirrors were maintained at liquid helium temper- 
atures inside the receiver to reduce radiative load- 
ing on the bolometers. The primary mirror, which 
was nutated to produce rapid azimuth modulation 
for MAXIMA, was fixed for MAXIPOL to avoid 
modulating the polarization properties of the tele- 
scope. 

Millimeter-wave radiation from the sky was re- 
imaged to a 4 X 4 array of photometers at the fo- 
cal plane. Observations were made in bands cen- 
tered on 140 and 420 GHz with Ai/ ~ 30 GHz. 
The twelve 140 GHz photometers were optimized 
to measure the CMB and the four 420 GHz pho- 
tometers were used to monitor foreground dust 
contamination. The 10' FWHM Gaussian beam 
shape for the 140 GHz photometers was defined by 
a smooth- walled, single-moded conical horn and a 
cold Lyot stop. The 420 GHz photometers em- 
ployed multi-mode Winston horns. The bolome- 
ters were maintained at 100 mK by the combina- 
tion of an adiabatic demag netization refrigerator 
(|Hagman fc Richarddfliosf) and a 300 mK ^He re- 
frigerator. A photograph of the focal plane and 
a cross-sectional overview of the optical system, 
which was essentially unchanged from MAXIMA, 
are shown in Figure [T] 

Figure [5] shows a cross-section of the portion 
of the cold optics that was modified to convert 
MAXIMA to MAXIPOL. We used a 3.175 mm 
thick A-cut sapphire HWP. Reflections from the 
HWP were minimized by bonding 330 fim thick 
wafers of Herasil to each face of the sapphire with 
Eccobond 24, an unfilled, low viscosity epoxy that 
was used to achieve glue layers as thin as 13 /im. 
The HWP thickness was selected to minimize the 
fraction of emerging elliptically polarized intensity 
and thereby optimize the overall modulation ef- 
ficiency of the polarimeter. The calculated effi- 
ciency incorporated the effects of the finite spec- 
tral bandwidth of the photometers and the oblique 
incidence of rays on the HWP. 

Since the anti-reflection coating was not bire- 
fringent, the two incident polarization orientations 
had different coefficients of reflection. This differ- 
ential reflection gave rise to a HWP synchronous 
signal at a frequency of 2/o, where fo is the ro- 
tation frequency of the HWP. To minimize this 



effect, we calculated the anti-reflection coating 
thickness that would minimize the difference in re- 
flection coefficients given the spectral bandwidth 
of the 140 GHz photometers, the thickness of the 
Eccobond 24 layer and the angles of the rays. The 
2/o signal was not a source of systematic error be- 
cause it was out of the polarization signal band 
centered on 4/o (see Section [5]). 

The polarization analyzer was a commercial 
grid polarizer epoxied to a rigid "roof-shaped" 
frame that was positioned over the horn open- 
ings (see Figure [S]). This roof-shaped polarizer 
was implemented so that the light that was re- 
flected by the grid was directed out of the opti- 
cal path and into a cold millimetcr-wavc absorber 
()Bocklll994l) i.hat was mounted on both sides of 
the focal plane. The grid polarizer had 98 electro- 
formed gold stripes per cm that were bonded onto 
a 38 /im thick Mylar substrate. The stripes were 
5 /im wide. 




Focal Plane 

Fig. 3. — A sketch of the roof-shaped analyzer positioned 
in front of the array of horns at the focal plane. 

The HWP was located near the Lyot stop and 
was driven from its axis at fo = 1.86 Hz during 
all observations (see Figure [2|) . The orientation of 
the HWP was measured with a 17-bit optical en- 
coder. This rotation frequency, combined with the 
azimuth scan frequency that was either 0.10 Hz 
or 0.06 Hz during the primary CMB observations, 
gave approximately three to five modulations of 
the Q and U Stokes parameters per 10 arcmin sky 
beam in one azimuthal scan. 

3. Observations and Scans 

The instrument was launched at 15:14 UT (9:14 
AM Local Time) on 2003 May 24. The first ob- 
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servation began at 18:08 UT when the payload 
reached an altitude of 38.7 km. The flight ter- 
minated 26 hours after launch, 477 km west and 
120 km south of the launch site. In this paper 
we only discuss the data collected during one 7.6 
hour-long nighttime CMB observation near Beta 
Ursae Minoris (/3-UMi), RA = 14h 50' 42.5", Dec= 
-1-74° 09' 42." O ther observations are discussed in 
JohnsonI (|2004h andlCoUind (|2006h . 
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Fig. 4. — Map of the total number samples per 3' pixel 
from the twelve 140 GHz polarimeters for the /3-UMi ob- 
servation. The white square marks the 2.3° X 2.3° region 
used for power spectrum estimation. 



Jupiter was observed during both the day and 
nighttime portions of this flight to obtain the abso- 
lute intensity calibration of the instrument, to map 
the telescope beam shapes, and to measure the 
location of the beam centers for pointing recon- 
struction. During these observations the gondola 
was scanned in azimuth with a triangie-wave pro- 
file centered on Jupiter. The turnarounds on this 
profile were smoothed to ensure quiescent noise 
performance from the bolometers. The daytime 
scans were ±0.6° with a 20 sec period. The night- 
time scans were ±1.3° with an 18 sec period. Si- 
multaneously, the telescope elevation was rastcred 
±0.5° relative to the planet every 20 min. 

For the CMB observations the telescope az- 
imuth was scanned in a similar fashion near (3- 
UMi, but with an amplitude of 1° and a period 
that was changed mid-scan from 10 to 15 sec. To 
improve cross-linking the telescope elevation was 



re-pointed 0.2° above and below the guide star ev- 
ery 10 min. Figure |4] shows the total number of 
140 GHz detector samples per 3 arcmin pixel pro- 
duced by this scan strategy. Observation statistics 
including pointing error for all scans are given in 
Table [H 

4. Attitude Reconstruction 

Telescope pointing was primarily determined 
with one of two star cameras. The camera used 
during daytime observations was filtered with a 
695 nm Schott glass filter to reduce the brightness 
of the daytime sky, and fitted with a reflective lens 
that had a 500 mm focal length, which provided 
a FOV of 0.72° by 0.55°. The unfiftered night- 
time camera used a 50 mm lens that provided a 
7.17° by 5.50° FOV. Pointing reconstruction us- 
ing star camera d ata has been described in previ- 



^pr 

ous publications ( Hanany et al. l l2000HRabiill2002 : 
20061: Icollinsll2006^ . 



Rabii et al 



MAXIPOL adopts the convention for the 
Stokes parameters /, Q, and U used by WMAP, 
which references the polarization direction to the 
NOP (jHinshaw et all l2003l : IZaldarriaga fc Seliak 
19971 ). When the instrument observes pixel p at 



time t, the output of the detectors (in units of volts 
and ignoring noise and systematic errors that wifl 
be discussed later) can be written 



dplt) = R[Ip 



e Qp cos ( 2'y{t) ) 
e Up sin ( 2-f{t) ) ] 



(1) 



Here 7 is the angle between the polarization ref- 
erence vector at sky pixel p and the polarimeter 
transmission axis, e is the polarimeter modulation 
efficiency, and R is an overall calibration factor 
that has units of V K"! (|Collins.,2006i) . The angle 
7 is given by 



7 a - 2/3- 



(2) 



Here, a is the angle between the reference vec- 
tor at pixel p and a vector pointing from p to the 
zenith along a great circle, (3 is the HWP orienta- 
tion angle in the instrument frame, and the 7r/2 
accounts for the fact that the transmission axis of 
the analyzer was perpendicular to the zenith axis. 
The angle a was computed for each time sample 
from the pointing reconstruction. The angle /? was 
computed by subtracting a photometer-dependent 
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Pointing Error 

Scan Target Start Time Duration Random Systematic 
[UT] [hours] [arcmin] [arcmin] 
Daytime Jupiter 20:03 L02 (199 Ol 
Nighttime Jupiter 02:42 0.52 0.47 0.21 
/3-UMi 03:18 7.64 0.32 0.27 

Table 1: Scan targets with scan durations and pointing error. A detailed description of the pointing errors is given in lCoIIin3 
l|2006r i. 



offset from the HWP encoder angle. This encoder 
offset was measured during laboratory calibration 
(see Section 15^ . By combining Equations[T]and[2] 
the model for noiseless time ordered data (TOD) 
can be written 

d{t) =R[Ip - eQp cos( 4/3(t) - 2a{t) ) 

+ eJ/p sin( 4/3(0 -2a(t))]. (3) 

5. Time Ordered Data Processing 

In this paper we report on the data from the 
twelve 140 GHz polarinicters. The raw data 
from /3-UMi consists of 5.76 million time-ordered 
samples for each of the 140 GHz polarimeters. 
The bolometer sample period was 4.8 msec. We 
first give an overview of the time-domain data- 
processing algorithm, which is outlined in Fig- 
ure El and then discuss some of the processing 
steps in more detail. The rest of the data analy- 
sis procedure is described in the companion paper, 
Wu et all (l2006h . 

Transients such as cosmic ray hits and inter- 
nal calibration pulse s were flagged using an al- 
gorithm described in IJohnsonI (j2004[ ). The data 
were calibrated and blocks of data separated by 
gaps longer than 30 sec were processed as sep- 
arate segments. These time domain data con- 
tained a significant instrumental signal that was 
synchronous with the rotation of the HWP. We 
refer to this instrumental signal as the HWP syn- 
chronous signal (or sometimes in short the Q). An 
initial estimate of the fi was subtracted. Gaps 
shorter than 30 sec were filled with constrained 
noise realizations, and the f5 was replaced, leaving 
continuous, transient-free data. The instrumen- 
tal filters were deconvolved after padding the ends 
of each segment with 100 msec of matched white 
noise, simulated HWP synchronous signal, and a 
window function that smoothly decayed to zero 
at the endpoints. A section of data one filter- 



time-width long was removed from the ends of 
each data segment after filter deconvolution to 
avoid any contamination from edge effects pro- 
duced by the Fourier transforms. The HWP syn- 
chronous signal was then re-estimated and sub- 
tracted from the raw data, producing the time- 
ordered data (TOD). The time-ordered polariza- 
tion data (TOPD) were extracted from the TOD 
by demodulation using a phase-locked, sine-wave 
reference signal constructed from the angle 27. 
The TOPD were checked for noise stationarity us- 
ing a frequency-domain test and for Gaussian- 
ity using the Kolmogorov-Smirnoff (KS) test in the 
time domain. 
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Fig. 5. — A flowchart of the time-domain processing al- 
gorithm. Blue boxes at left are data inputs, numbered red 
boxes in the middle are analysis operations and green boxes 
at right are data outputs. 

The following subsections describe the details of 
the HWP synchronous signal estimation, the filter 
deconvolution, the TOPD demodulation, and the 
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tests of the noise properties. 

5.1. Subtraction of the HWP Synchronous 
Signal 

The raw data were dominated by the HWP 
synchronous signal, which came from thermal 
emission from the HWP drivetrain, differential 
transmittance from the HWP, and modulated 
instrumental polarization. Emission from the 
HWP drivetrain appeared at harmonics of fo with 
power which decreased with frequency. Differen- 
tial transmittance contributed primarily at 2/o, 
and polarized signals generated by instrumental 
polarization were dominant at 4/o. 

We modeled the fi as a sum of n = 8 harmonics 
of the HWP rotation frequency with amplitudes 
^„ and phases Xn 

8 

HWPSS {t) = ^" (*) sin ( "/^(t) + Xn (t) ) (4) 

n=l 

where HWPSS is the fi. Recall that (3{t) is the 
HWP rotation angle in the instrument frame. 
We used only eight harmonics because signals 
at higher frequencies were blocked by the fil- 
ters discussed in Section 15.21 Each of the three 
sources listed above (thermal emission, differen- 
tial transmission, and instrumental polarization) 
contributed its own amplitude and phase to pro- 
duce the net amplitude ^ and phase x each of 
the harmonics. We allowed the amplitude of each 
of the three sources to change linearly with time 
but we fixed their phase angles relative to the in- 
strument frame. Thus the contribution of the nth 
harmonic of the 6 can be written as 

HWPSS„ (t) = e„ (0 sin ( n/3(t) + Xn {t) ) 

3 

= ^ (A,y + Bnjt) sin ( nf3{t) + 0„j ) (5) 

where Anj, Bnj, and 4>nj are constants. Trigono- 
metric identities can be used to relate Equa- 
tions [4] and [5] and to rewrite the fi as 

8 

HWPSS (i) = ^ ( Ci„ + C2nt ) cos ( np{t) ) 

n=l 

+ {Sin + S2nt) sin (71/3(0) (6) 

with appropriate relations between ^, x. A, B and 
the constants Ci„, C2n, Sin, and S2n- We sub- 
tracted the HWP synchronous signal from the raw 



time domain data using fits to Equation [6] in the 
following way. 

The 32 C and S coefficients were assumed to 
be constant within each data segment and were 
estimated in each segment by an iterative demod- 
ulation procedure. To begin the first iteration, the 
data segment was demodulated for n = 1 with the 
reference signal cos(/3). The coefficients Cn and 
C21 were estimated using a linear least-squares 
fit to the output, ignoring the flagged data con- 
taminated with transients. The component cor- 
responding to Cii and C21 was constructed and 
subtracted from the data. The subtraction was re- 
peated with the sin(/3) reference signal, and then 
for the harmonics n = 2, ...,8 in order. The pro- 
cess was then iterated 25 times. After each itera- 
tion the new estimates of C and S were added to 
the old estimates. For a typical photometer, after 
four iterations of this process the difference be- 
tween the data and the HWP synchronous signal 
was observed to be indistinguishable from random 
noise (see Section lS.f .ip . 

The amplitude values ranged from approxi- 
mately 1.5 to 106 mK for harmonic n = 1, from 30 
to 250 mK for n =2, and from 33 to 600 mK for 
n = 4. The amplitude drifts were typically 0.5% 
over a 10 min data segment and 10% over a 3 hour 
time scale. The phases Xn typically varied by less 
than 5 deg over the entire CMB scan. More de- 
tails about the characteriza tion of the HW P syn- 
chronous signal are given in ICollind (I2OO6I) . 

Figure [6] shows the properties of one segment 
of data from one of the polarimeters before and 
after subtraction of the 6. Panel (c) shows that 
the power spectrum of the data after the fi was 
subtracted is flat in the vicinity of 4/o = 7.4 Hz 
and also down to frequencies well below 1 Hz. This 
white-noise level is the nominal noise level of the 
instrument. 

The subtraction of the fi was tested using sev- 
eral figures of merit which address the following 
questions. Are the data after subtraction Gaus- 
sian distributed? Does the subtraction of the fi re- 
move CMB signal from the map? Are there resid- 
uals of the fi in the map domain? How stable is 
the 4/0 component of the fi and is it well-fit by 
the slow time variation of the model? We will ad- 
dress the first three questions here and the others 
in Section [531 
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Fig. 6. — (a) A 16 mill long sample of raw data from one 140 GHz polarimetcr plotted versus the HWP angle /3. (b) Power 
spectrum of the raw data, (c) Power spectrum of the data after subtraction of the HWP synchronous signal (TOD), (d) Power 
spectrum of the demodulated TOD, which is the TOPD. Zero- frequency in (d) is 4/o in (b) and (c). The red, dashed curves 
in panels (fe), (c) and (d) mark the polarization signal bands limited by the telescope scan frequency and the beam. The blue, 
dash dot curves in panels (fe) and (c) mark the temperature signal band. Some 1// noise appears in the temperature band. 



5.1.1. Gaussianity 

The properties of the data after the HWP syn- 
chronous signal was subtracted were tested by 
band-pass fihcring the TOD between 0.3 and 8.9 Hz 
to select the band that contains the first four 
Fourier modes of the fi. The data were decor- 
related by filling the rejection band with a white 
noise realization, and then averaged into 36 bins in 
the HWP angle domain, each spaced by 10°. The 
10° bin size was selected so any residual 4/o signal 
could be resolved. A statistic was calculated 
over the 36 bins for each data segment to search 
for residual HWP synchronous signal. The set of 

from each polarimetcr was compared to the 
distribution for 36 independent degrees of freedom 
using the KS test, and no excess was found. 



5.1.2. Removal of CMB Signal 

We simulated the removal of the HWP syn- 
chronous signal with simulated raw time streams 
that had the same properties as our real time 
streams. The fi, the pointing information, and 
the noise estimate from each segment of data were 
used with simulated CMB polarization maps to 
assemble a simulated raw data stream composed 
of the fi, the HWP-modulated CMB polarization 
signal, and a noise realization with a power spec- 
trum given by Equation [9l The subtraction of the 
fi was performed on each segment, the noise was 
subtracted, and the CMB signal was recovered by 
demodulation (see Section [5T3|) . The initial CMB 
TOPD was subtracted from the final. The residu- 
als were bandpass filtered between 0.05 and 1.5 Hz, 
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and then averaged into maps of 3' pixelization us- 
ing weights set by the noise level in each segment. 
In the 2.3° x 2.3° square region used in the power 
spectral analysis the residual maps had an RMS 
of 0.018 ± 0.002/xK in Q and 0.016 ± O.OOl^iK in 
U over five trials. These residuals were negligible 
compared to the RMS of simulated CMB maps of 
the same region, which was 4 /iK for both Q and 
U. 

5.2. Filter Deconvolution 

The TOD were filtered in the readout electron- 
ics with high- and low-pass Butterworth filters 
with edges at approximately 15 niHz and 20 Hz, 
respectively. In addition, the frequency response 
of the bolometers attenuated and phase-delayed 
the modulated 4/o signals at 7.4 Hz by about 
35% and 20° of HWP angle, respectively. The 
frequency response of the bolometers was mea- 
sured in the laboratory. The results were modeled 
with two time constants and the parameters are 
shown in Table [21 Alternative analysis methods 
gave somewhat different parameters for the values 
of the time constants and of the relative weights. 
These differences were larger than the errors on 
the parameters within any given method. How- 
ever, calculations showed that the differences be- 
tween any of the derived filter functions in terms of 
amplitude and phase response at the signal band- 
width near 7.4 Hz were small such that the ef- 
fect on the final results would have been negligible 
compared to the calibration uncertainty or to the 
uncertainty on mixing between the E and B modes 
due to noise. 

The effects of the electronic filters and bolome- 
ter time constants were removed from the TOD 
by deconvolving the complex filter J^, where 

^{f) = ^hp^lpJ^bolo/\^lp3o\- (7) 

Here, J-hp and J-ip arc the electronic high and low- 
pass filters, J-boio is the bolometer response func- 
tion, and |.77p3o| is a real, phase-preserving soft- 
ware filter with a low-pass cutoff of 30 Hz. This 
filter was included to suppress the spurious high 
frequency noise created by the deconvolution of 
the bolometer response function. 

5.3. Demodulation 

Modulation by the rotating HWP, combined 
with telescope scanning, moves the polarized sky 
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Fig. 7. — Top panel: The power spectrum of a 12.4 min 
long segment of TOPD for Q with the binned power spec- 
tra of the first and second halves overplotted (red stars and 
blue diamonds) . The binned power spectra are slightly off- 
set in frequency from the bin center for clarity. This seg- 
ment was deemed stationary by the stationarity test, with 
= 12.52 for 9 degrees of freedom. The cutoff ^™ 
stationarity for the 12.4 min long chunk length was 19.43 
(see Section l5.4l l. Bottom panel: the difference between the 
cumulative distribution function (CDF) of decorrelated Q 
data from one data segment and the CDF of a Gaussian dis- 
tribution of the same standard deviation. The KS statistic 
showed no non-Gaussianity above the 0.26(t level for this 
segment. The red dashed lines mark the Icr significance 
level. 

signals to sidebands of 4/o, as shown in Figure [6l 
To extract the Q and U sky signals we used a soft- 
ware lock-in technique. The TOD were demodu- 
lated with the reference signals — cos(4/3 — 2a) /e 
and sin(4/3-2a) /e to produce the Q and U TOPD, 
respectively. Demodulated polarized sky signals 
lie in a band defined at the low frequency side by 
the '--^O.OB Hz telescope scan frequency and at the 
high side by attenuation from the beam above ap- 
proximately 1.5 Hz, see Figure[6l For map making, 
the demodulation procedure also included a band- 
pass filter with high and low-pass edges at 0.05 
and 1.5 Hz, respectively, that was used to reject 
any out-of-band signals. 

Figure [Gjd) shows the power spectrum of one 
segment of Q TOPD. This power spectrum is flat 
down to a frequency of approximately 1 mHz. The 
flatness is a measure of the degree of stability of 
the 4/o component of the HWP synchronous signal 
and the efficacy of the subtraction method. For 
99.5% of the TOPD from all photometers the 1// 
knee was below 0.06 Hz. 
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Photometer 


NEQ,U 


Modulation 




Wa, Wb 


Calibration 


FWHM 






Efficiency 


[ms] 


[weight] 


xlO-5 [V K-i] 


[arcmin] 


bl3 


780 


0.92 


5, 27 


0.78, 0.22 


2.30 


9.2, 10.2 


bl4 


2000 


0.93 


11, 41 


0.55, 0.45 


1.56 


9.8, 9.9 


bl5 


1100 


0.90 


8, 29 


0.64, 0.36 


2.32 


9.3, 10.8 


b23 


780 


0.93 


8, 41 


0.66, 0.34 


2.18 


10.0, 11.1 


b24 


870 


0.93 


8, 50 


0.28, 0.72 


3.32 


10.2, 11.3 


b25 


710 


0.93 


10, 56 


0.44, 0.56 


2.91 


10.3, 10.7 


b33 


1000 


0.93 


9, 46 


0.50, 0.50 


2.15 


9.4, 10.4 


b34 


910 


0.93 


8, 26 


0.39, 0.61 


1.98 


9.3, 10.6 


b35 


880 


0.94 


12, 64 


0.49, 0.51 


2.34 


9.5, 10.6 


b43 


1200 


0.95 


9, 46 


0.50, 0.50 


1.93 


8.6, 11.3 


b44 


1700 


0.93 


11, 74 


0.38, 0.62 


1.90 


9.0, 11.4 


b45 


1100 


0.92 


14, 57 


0.42, 0.58 


2.42 


9.8, 12.1 


receiver 


280 


0.93 


9, 46 


0.50, 0.50 


2.28 


9.5, 10.9 



Table 2: Characteristics of each of the 140 GHz polarimeters in the array. The row labeled "receiver" contains the inverse square 
sum of all NEQ,U values thereby giving the total receiver performance. Other elements in this row are simple averages. The 
two bolometer time constants, Tj\ and Tg, and their corresponding weights, Wa and Wbi are differentiated with the subscripts 
A and B. Typical beam FWHM uncertainties from the MCMC analysis are 2%, while typical calibration and modulation 
efficiency uncertainties are 13% and 2%, respectively. The NEQ,U and the calibration are consistent with an instrument 
model that includes the bolometer noise, the amplifier noise, the measured bolometer time constants, the transmission of the 
analyzer, and the reflectance of the HWP. 



5.4. Noise Characterization 

The TOPD were checked for stationarity by 
bisecting each data segment and comparing the 
power spectra of the two halves. For this com- 
parison, the portion of the power spectra that 
spanned the polarization signal band between 
0.02 and 1.5 Hz was divided into nine bins. This 
choice of binning gave at least 8 frequency modes 
per bin for the shortest data segments, which were 
approximately 2 min long. The following statis- 
tic 

was used to assess the similarity of the power spec- 
tra. Here, the sum extends over the nine bins, 
and Pi and Ui were the mean and standard er- 
ror of the power spectrum in each bin, respec- 
tively. The distribution of this statistic for white 
noise time streams was estimated with Monte- 
Carlo simulations of 20,000 realizations. Seven dif- 
ferent segment lengths ranging from 2 to 24 min 
were considered. The range of segment lengths 
were selected based on the maximum and mini- 
mum lengths of the data segments used for map 
making. The corresponding to the probabil- 



ity to exceed 2a was stable to 3% over this range 
of segment lengths and was estimated for an ar- 
bitrary segment length by spline-interpolating the 
2cr between the seven simulated lengths. 

The criterion of 2a stationarity was used with 
an algorithm that searched for long sub-segments 
of stationary data. A routine recursively cut each 
data segment into two halves if the two power 
spectra were dissimilar. After bisection, an it- 
erative routine attempted to concatenate adja- 
cent sub-segments whose power spectra were not 
yet known to be dissimilar. The sub-segment 
was flagged as non-stationary if the sub-segment 
length after this procedure was less than 2 min. 
Non-stationarity was mostly confined to three of 
the twelve 140 GHz photometers for which approx- 
imately 11% of the TOPD were discarded as non- 
stationary. The average stationary sub-segment 
length was 9 min. The upper panel of Figure [7] 
shows an example of a segment of data that was 
found to be stationary. 

To quantify any 1// noise in the TOPD the 
following three parameter model 



P{f) = 



1 + 


/ fknee \ 




\ f ) \ 



(9) 
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was fit by least-squares to the binned Pi of the 
stationary segments with weights ai. The model 
parameter fknee is the 1// knee frequency and a is 
the spectral index of the low-frequency noise. The 
white-noise parameter tr^ was used to calculate 
the noise equivalent values of Q and U shown in 
Table [2] by averaging over all stationary segments. 

The decorrelated time streams Qd and Ua were 
produced by band-pass filtering the TOPD be- 
tween 0.02 and 1.5 Hz and replacing the out of 
band frequencies with white noise realizations in 
the frequency domain. The KS test was performed 
comparing the cumulative distribution function 
(CDF) for Qd and Ud to the CDF of a Gaus- 
sian distribution of standard deviation equal to the 
standard deviation of each segment. All station- 
ary segments passed the KS test to better than 2a 
confidence, where the confidence level is referred 
to the normal distribution. The bottom panel of 
Figure [7] shows an example of this procedure for 
one segment of the data. 

6. Calibration 
6.1. Responsivity 

MAXIPOL was calibrated with observations of 
Jupiter, which provided the beam shapes, the pho- 
tometer responsivities, and the beam centers for 
pointing reconstruction. In addition, the Q and 
U beams were mapped to verify the polarization 
properties of the instrument. The time ordered 
data had contributions from the planet signal, 
transients, the HWP synchronous signal and noise 
with a 1// component. We used the procedure de- 
tailed in the Appendix to make a binned map of 
Jupiter for both the day and nighttime scans for 
each of the photometers. 

The beams were characterized by fitting a 
Gaussian to the data with a Bayesian likelihood 
analysis that used a Monte Carlo Markov-Chain 
to explore the likelihood function. We used this 
method to mitigate the effect of unsampled pixels 
in the maps, as the method explicitly considers 
errors in each pixel. We determined the follow- 
ing beam parameters: the overall amplitude, the 
width of the beam along two axes, the orienta- 
tion angle of the ellipticity, and an overall offset, 
as well as the location of the photometer in the 
focal plane. In addition to determining beam pa- 
rameters from the binned map, we also applied 



the analysis directly on the time ordered Jupiter 
data. The results for the two methods were consis- 
tent within the errors. The resulting beam widths 
are reported in Table [21 and the beam models are 
plotted in Figure [S] The asymmetry of the beams 



was acc ounted for using the recipes in IWu et al 
(|2001b[ ). 



The calibration R reported for each photome- 
ter in Table [2] was computed by integrating the 
maximum-likelihood Gaussian beam 



R 



1 



V{x,y) dx dy (10) 



where 



S'{,^) (dB^j/dT) dv 

T=^ . (11) 

j S{u) (dBcMB/dT) dv 

Here, Aj is the solid angle of Jupiter, i?cMB is 
the Planck function, Sj^j is the Rayleigh-Jeans 
brightness, V{x^y) is the beam map in volts and 
S and S" are the spectral response of MAXIPOL 
and the calibration apparatus, respectively. The 
Rayleigh-Jeans brightness temperat ure of Jupiter 
was taken to be Tj ~ 173 ± 9 K (jCriffin et al. 
Il986ll . 

The temperature of the detector assembly was 
not regulated so the temperature of the detec- 
tors rose with time at a rate of approximately 
7 mK hr~^. In order to monitor the time depen- 
dence of the calibration, the bolometers were illu- 
minated by a fixed-intensity millimeter- wave lamp 
for 10 sec every 22 min. The resulting bolome- 
ter signals were found to decrease linearly with 
increasing bath temperature. The relative respon- 
sivities measured during the flashes of the lamp 
were used to extend the absolute Jupiter calibra- 
tion to all times during the observations using the 
measured temperature of the detector assembly. 

6.2. Polarimeter Characterization 

We define the modulation efhcicncy e as the 
multiplicative factor that relates the measured 
data to the incident Q and U Stokes compo- 
nents, see Equation [3l It can be shown that for 
monochromatic light of frequency v 



(.y = g sm 



(12) 
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Fig. 8. — (a) A beam map for the photometer bl4 from the daytime observation of Jupiter, (b) Contour Unes from the 
maximum hkehhood Monte Carlo Markov-Chain (MCMC) fits for all the 140 GHz photometers. Panel (a) shows the 1, 10, 50 
and 90% contour lines. The smooth lines come from the Gaussian fit and the jagged lines from the raw data. Panel (b) shows 
the 10, 30, 50, 70 and 90% contours from both MCMC analyses discussed in Section |6. II 



where (S(i^) is the retardance of the HWP and g is 
the polarization efficiency of the wire-grid polar- 
izer (|johnsonll20"o3 : ICollinsll2006f) . We calculated 
a typical expected efficiency 



tvS{v)dv 
S{v)dv 



(13) 



from the thickness of the HWP, the measured g, 
and a typical spectral response as a function of fre- 
quency S{v). For MAX IPOL g was measured to 
be 0.97 (|johnsonl[200l . The expected efficiency 



Ceajp was 93%, and this result agrees with the effi- 
ciencies measured in the laboratory and listed in 
Table[2]to within the 2% experimental uncertainty. 
The following section explains the laboratory mea- 
surement and analysis. 

A wire-grid polarizer made from the same ma- 
terial used to make the polarization analyzer was 
mounted on the cryostat window with its trans- 
mission axis oriented to within 2° of the axis of 
the analyzer. Beam filling thermal radiation from 
a 273 K ice bath was chopped at approximately 



6.5 Hz with a 300 K aluminmn chopper blade 
covered with absorbing foam (0.64 cm thick Ec- 
cosorb LS-14). In order to avoid bolometer sat- 
uration from the brightness of the warm loads, 
a 4 K absorptive attenuator (1.9 cm thick Ec- 
cosorb MFllO) that was designed to transmit less 
than 3% at 140 GHz was inserted into the op- 
tical p ath at the intermediate focus of the tele - 
scope (jPeterson fc Richardsl[l983 |johnsonll2004[) . 
The HWP was stepped in ~ 5° intervals and ap- 
proximately 20 sec of data were collected at each 
HWP orientation. Detector drifts were removed 
by fitting a model consisting of a second-degree 
polynomial and a single-timc-constant exponen- 
tial to each data segment. The transmitted chop 
amplitude was estimated by subtracting the noise 
RMS, which was obtained from an independent 
noise-only measurement, from the data RMS. The 
effect of bolometer bath temperature drifts was 
mitigated by implementing a time-dependent re- 
sponsivity correction. The detector response was 
linearized by applying a quadratic, amplitude de- 
pendent responsivity correction. 

The data from one typical photometer are 
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Fig. 9. — Data from a laboratory measurement of the 
modulation efficiency. The setup for this measurement is 
described in Section 16.21 The solid curve is the best fit 
model given by EQuation ll4l The reduced for this model 
is 1.03 for 42 degrees of freedom. From the fit parameters, 
we calculated V to be 0.88 ±0.03, which corresponds to po- 
larimeter efficiency e = 0.94 ±0.02. The difference between 
these data and the best-fit model are plotted in the bot- 
tom panel to show the goodness-of-fit and the measurement 



shown in Figure [9l The foUowing five parame- 
ter model was fit to the binned data 



m(p) = A2/„ cos(2/3 + 02/„ 
+ A^f^ cos(4/3 + 04 



B. 



(14) 



The nominal chop amplitude error per HWP ori- 
entation bin was estimated to be the detector noise 
RMS. This random error was assigned before re- 
sponsivity corrections. The subsequent responsiv- 
ity corrections caused the magnitude of the ran- 
dom error to vary from bin to bin. A system- 
atic error was estimated by subtracting the best 
fit model from the binned data, computing the 
RMS of the residual and then adding this error in 
quadrature with the nominal error per bin. This 
combined error was used to compute the final er- 
rors in the fit parameters. The experimental figure 
of merit 



(15) 



was calculated from the fit parameters using the 
equivalent expression V = A4f^/B, and the error 
in V was propagated from the errors in ^4/„ and 
B. Using Mueller calculus it can be shown that V 



l + .g2 V 
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(16) 



and the value for e for each photometer is given 
in Table O The agreement between e and the ex- 
pected modulation efhciency eexp (Equation I13p 
confirms that the 7% average depolarization for 
the 140 GHz photometers is accounted for by the 
efficiency of the grid polarizer and the spectral re- 
sponse of the HWP. 

7. Analysis and Results 

We used a maximum-likelihood method to 
make maps of Q and U from the TOPD (see 
Figure \W\\ . The map pixel size was set to 3 ar- 
cmin. Two approaches, one frequentist and one 
Bayesian, were used to estimate the EE, BB and 
EB angular power spectra in three £ bins: £ < 150, 
151 < ^ < 693 and £ > 694. Given the beam size 
and the sky coverage we expect that only the cen- 
ter bin would have any signal. For both methods 
we used the 2.3° x 2.3° region of the maps inside 
the white square in Figure [H The average integra- 
tion time in this box was 117 sec for each 3 arcmin 
pixel. A co mprehensive des cription of the analysis 
is given in IWu et al.l ( 2006f ). Here we only give a 
summary. 

For the Bayesian analysis we estimated the 
amplitude of the angular power spectra us- 
i ng a Monte-Carlo Markov Chain algorithm 



(jMacKavl l2003l : iLewis fc Bridie l2002h derived 



from the MADCA P spectrum solver MADspec 
( Borrill et al ]l2006l) . -^Tom a pair of 50,000 element 
chains we calculated the posterior likelihoods of 
the amplitudes of the power spectra in all 9 £ bins. 
Estimates for the amplitude were calculated both 
with and without priors which set the EB and BB 
spectra to zero with four different functions that 
describe the shape of the power spectrum within 
the bin. Here we quote the results for the ampli- 
tude of £{£ + l)Ci/2TT after marginalization over 
the eight un-interesting £ bins and without any 
priors on the EB and BB power spectra. Results 
obtained with o ther assumptions are presented in 
(120061 ). 



Wu et al 



The left panel of Figure [TT] shows the esti- 
mated maximum-likelihood EE power and other 
reported results, while the right panel shows the 



13 



111 113 115 111 113 115 111 113 115 111 113 115 

Z (degrees) /(degrees) /(degrees) /(degrees) 

^^^^^M : : I ^iK ^^^^^fc : i 

-200 -100 100 200 -20 20 



Fig. 10. — Q and U maps of the MAXIPOL observation region near /3-UMi created from the raw data (left) and raw maps 
convolved with a 10 arcmin FWHM Gaussian beam (right). The pixel size is 3 arcmin. The convolved maps are used for visual 
inspection only. 



corresponding posterior likelihood function. This 
result gives weak evidence for an EE signal with 
a maximum-likelihood amplitude of 55I45 t^^^ 
(68%). The likelihood function is asymmetric and 
skewed positive such that if we assume a uniform 
prior over all possible values of £{£+l)Ci/2TT (both 
positive and negative) the probability for a signal 
larger than zero is 96%. Including the calibra- 
tion uncertainty the maximum-likelihood ampli- 
tude is 53I45 /^-^^ (68%) with no change in the 
probability of positive power. The EE amplitude 
pr edicted by the concordance ACDM model given 
bv ISpergel et all (|2006h is 14 ^K^ This theoret- 
ical value falls inside the 65% confidence region 
of our likelihood function around the maximum- 
likelihood value. The sharp cut-off on the nega- 
tive end of the parameter axis is a consequence 
of correlations that occur within the high dimen- 
sionality of the marginalized space. Despite this 
feature, the likelihood accur ately represents the 
result of the experiment (see Wu et al. I 2006I ) for 
more details). 

Both the 68% confidence intervals and the sig- 
nificance of positive power depend somewhat on 
the shape function used during the power spec- 
trum estimation and on whether there are prior 
constraints on the EB and BB spectra. For ex- 
ample, assuming a uniform prior, the probability 
that the EE power spectrum amplitude is positive 
is 83% for a shape function of !/£{£+ I) with EB 
and BB set to zero. The probability is 98% for 
a shape function of l/{2£+ 1) with no constraints 
on EE and EB. 



The amplitudes of the EB and BB power spec- 
tra assuming the l/i{£ + 1) shape function are 
IStll and -3lt?9 ^iK^ (68%), respectively, which 
are consistent with ACDM predictions. MAX- 
IPOL does not have the sensitivity to detect either 
BB or EB. 

For the frequentist analysis we computed 
pseudo-band powers from the Fourier transform of 
the Q and U maps using the flat sky approxima- 
tion. The pseudo-band powers have biases coming 
from the beam convolution, the time domain pro- 
cessing, and the Fourier transform. These biases 
were corrected, and band-power error bars were 
estimated using Monte Carlo simulations and the 
measured beam profiles. Results from the fre- 
quentist methods are consistent with those from 
the Bayesian analysis although the error computed 
by the Bayesian method is smaller, as expected. 

8. Systematic Errors 
8.1. Maps and Spectra 

We tested the Q and U maps used in the 
power spectrum analysis for signal Gaussianity by 
analyzing the eig envalue- normalize d Karhunen- 
Loeve coefficients (|Wu et al.l l2Q01al) . If the sig- 
nal is Gaussian then these coefficients should 
be normally distributed. Some of the eigen- 
values of the noise-whitened signal matrix were 
negative because of the high noise and imper- 
fectly estimated signal in those modes. We ex- 
cluded these modes from the test, but included 
noise-dominated modes. The resulting coefficients 
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Fig. 11. — The MAXIPOL estimate of EE power and 
other reported results (see Section [T] for references). The 
solid black line is th e ACDM concordance model given by 
ISpergel et aD l|2006l ). The MAXIPOL result shown here 
is the 68% probability region around the maximum of the 
likelihood function without any priors on the EB or BB 
spectra. The likelihood function is shown in the right panel. 
The black arrow in this likelihood panel marks the expected 
ACDM signal for our bin, which is 14 /xK^. Assuming a 
uniform prior the probability that the amplitude is posi- 
tive is 96%. Including the calibration uncertainty does not 
change the significance of positive power and increases the 
error bars slightly, see text. 

passed the Kolmogorov test for normality at 95% 
confidence. 

To check for residual systematic errors we dif- 
ferenced maps from the first and second halves 
of the TOPD. The power spectra calculated from 
these difference maps give a result that is consis- 
tent with the absence of signal and with random 
noise. The maximum-likelihood value for the EE 
power is -19tlt^lK^ (68%). 

8.2. Cross-Polarization 

We refer to "cross-polarization" as a leakage be- 
tween Q and U states of polarization as quantified 
by the QU cross term of the appropriate Mueller 
matrix. Here, we describe characterization of the 
cross-polarization of the main beam. 

Laboratory measurements outlined in Sec- 
tion 16.2! characterized effects due to the photome- 
ter dependent cross-polarization of the optical sys- 
tem, but without the effects of the primary mirror, 
combined with a photometer independent HWP 
encoder offset. Both effects were corrected in the 



analysis by using the photometer-dependent HWP 
angle (3. 

For the purpose of instrument characteriza- 
tion these two effects were separated by assuming 
that the center of the focal plane has zero cross- 
polarization because it is in the symmetry plane 
of the optical system. This assumption was sup- 
ported by simulations of the optical system. The 
variation in the /3 correction around this zero-point 
was interpreted as photometer dependent cross- 
polarization. The measurements show that the 
cross-polarization produces less than 3.4° of lin- 
ear polarization rotation for any photometer. This 
result agrees with ray-tracing simulations. 

Simulations show that adding the primary mir- 
ror to the system should increase the linear polar- 
ization rotation induced by cross-polarization by 
less than 0.1°. This additional rotation was not 
corrected for during data analysis because calcu- 
lations showed the level of error in the power spec- 
trum produced by the effect was negligible. 

Uncertainties in the angles a and (3 (see Sec- 
tion |4|) are predominantly small systematic offsets 
that are equivalent to unknown cross-polarization. 
The uncertainty in a is 2° due to the uncertainty in 
the orientation of the balloon gondola during the 
pre-flight telescope alignment. The uncertainty in 
P is less than 2° due to the uncertainty in the 
alignment of the wire grid on the cryostat win- 
dow during the polarization calibration, and due 
to the uncertainty in the knowledge of the bolome- 
ter time constants. According to simulation, these 
errors have a negligible impact on the results. 

8.3. Instrumental Polarization 

We refer to "instrumental polarization" as 
those effects that would produce a detection of 
polarized light even if an unpolarized light is inci- 
dent on the instrument. With HWP polarimetry, 
polarized light that is produced on the sky side of 
the HWP is modulated and gives rise to a system- 
atic signal at 4/o. The systematic signal can leak 
into the signal bandwidth, which is at a sideband 
of 4/o, if it varies with time. 

In MAXIPOL several mechanisms contributed 
to an instrumental polarization signal at and near 
4/o. Differential reflection of light by the pri- 
mary, secondary and tertiary mirrors and by the 
polypropylene vacuum window of the cryostat pro- 
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duced partially polarized light in transmission. Of 
these, the reflection by the window was dominant. 
Emission from the mirrors and from the window 
produces partially polarized light because of the 
asymmetric nature of the optical system. Of these 
emission from the primary mirror was dominant. 
Unpolarized emission inside the cryostat that is 
diffracted around sharp edges into the light path 
can produce substantial polarization. 

The instrumental polarization produced by 
these mechanisms was divided into two types: a 
stable polarized offset and a time varying instru- 
mental signal that was synchronous with the sky. 
The stable offset appeared as the 4/o component 
of the HWP synchronous signal that we discussed 
in Section 15.11 The photometer-dependent mag- 
nitude of this signal was in the range of 30 to 
600 mK across the focal plane (see Section 15. ip . 
Calculations show that instrumental polarization 
produced in transmission by differential reflection 
off the bowed polypropylene window of the cryo- 
stat will produce signals at the low end of this 
range. The upper end of this range is consistent 
with signals produced by diffraction from sharp 
edges of apertures within the cold optics box. 
This offset was stable on time scales much longer 
than data segments that were used for map mak- 
ing and therefore amplitude drifts did not leak 
into our signal band between 0.06 and 1.5 Hz (see 
Section [5.11 and Figure [Hd)). This stable signal 
was rejected both during subtraction of the HWP 
synchronous signal and during demodulation. 

The time- varying instrumental signal was mod- 
ulated with the intensity pattern on the sky, 
thereby leaking / — > Q, J7 and therefore T ^ E, B. 
It arose only through polarization produced in 
transmission. We measured this signal in two 
ways. The level of instrumental polarization pro- 
duced by the receiver alone was determined to be 
less than approximately 1% for a typical 140 GHz 
photometer in the laboratory before flight. The 
procedure used for this measurement was essen- 
tially identical to the one outlined in Section [6?2l 
However, for this measurement the polarizer was 
removed from the cryostat window so the chopped 
signal observed an unpolarized load. During flight, 
the effect was measured more accurately with 
the full instrument during the Jupiter calibration 
scan. Since the level o f polarization of Ju piter 
is expected to be small ( Clemens et al. 1990l ) the 



planet served as an unpolarized point source. The 
instrumental polarization produced in transmis- 
sion was measured by rotating the HWP during 
the beam-mapping procedure. The resulting Q 
and U beam maps of Jupiter yielded no detectable 
instrumental polarization signal above about 1% 
for ten of the twelve 140 GHz photometers. The 
remaining two photometers, which were located 
at the edge of the focal plane, detected instru- 
mental polarization at the level of 4% and 5%. 
Calculations show that this level of instrumental 
polarization is plausible given a particular align- 
ment between the focal plane and the window to 
the cryostat. Assuming conservatively that all 
photometers had approximately 4% instrumen- 
tal polarization in transmission, simulations show 
that this performance would only produce a spuri- 
ous 3 /xK^ signal in our EE and BB power spectra 
for the multipole bin oi £ = 151 to 693. Since this 
leakage signal is undetectable in our data, no cor- 
rection was applied during the analysis. 

8.4. Foregrounds 

The region of the sky near /3-UMi was selected 
for CMB observations because contamination from 
dust, synchrotron and point sources is expected to 
be negligible. Polarized dust was our main con- 
cern because synchrotron radiation is brightest at 
frequencies below the MAXIPOL spectral bands, 
and point sources should not be polarized. Ex- 
tr apolating from 100 /xm using Finkbeiner model 



8 (Finkbeiner et al 



1999f l. the unpolarized mean 
dust brightness over the 2.3° x 2.3° square re- 
gion used for power spectrum estimation should 
be 4.1 fjK with an RMS of 0.8 ^K. Archeops 
found a 5% polarized fraction for dust e mission 
in the galactic plane (jBenoit et al ] |2004J l. The 
full sky WMAP 3-year data are consist ent with 
this polarized fraction ( Page et al. 2006f l. Given 



these measurements, the expected foreground po- 
larization anisotropy over the /?-UMi square region 
should be less than 0.04 /iK RMS, which is not de- 
tectable by MAXIPOL. A catalog search yielded 
no detectable radio o r infrared sources in our fleld 
(jSokasian et al.ll200l[ ). No foreground corrections 
were made during data analysis. 

The data from the 420 GHz polarimeters can be 
used to improve the knowledge of the level of fore- 
ground dust contamination in the /3-UMi region. 
However, this 420 GHz information was not used 
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in this analysis for two primary reasons. First, 
the 420 GHz maps of Jupiter were poorly sampled 
so more work is required to determine the inten- 
sity calibration and the beam eenters for pointing 
reconstruction. Second, the attenuator used for 
pre-flight polarimeter characterization (see Sec- 
tion 16. 2p did not transmit detectable amounts of 
420 GHz signal, so a dedicated 420 GHz post- 
flight polarimeter characterization is required to 
calibrate the TOPD of this data set. 

9. Discussion 

MAXIPOL was designed to be a pathfinder for 
bolometric CMB polarization experiments. It is 
the first bolometric CMB experiment to report re- 
sults with a rapid polarization modulator. 

The predicted benefits of the HWP technique 
have been demonstrated. Power spectra of the 
time domain data from the twelve independent 
polarimeters gave white noise after demodulation 
to frequencies well below 50 mHz for most of the 
data. A significant fraction of demodulated data 
had white noise at frequencies as low as 1 mHz. 
The Q and U data showed no detectable system- 
atic errors despite a sizable HWP synchronous sig- 
nal in the raw data. There was no need to com- 
pare noise and responsivity data between detec- 
tors, which is necessary when using differencing 
polarimeters. This made the analysis simpler. 

We used the Q and U data to construct a map 
and estimate the polarization power spectra. The 
time domain noise was shown to be Gaussian and 
stationary. The noise in the map is consistent with 
Gaussian random noise at a level expected given 
instrument noise. The time streams and maps 
were subjected to multiple tests for systematic er- 
rors with null results. 

The data give weak evidence for EE power at 
a level that is consistent with the prevailing cos- 
mological model. The EB and BB power spectra 
are consistent with zero, also as expected given 
instrument noise and the cosmological model. 

MAXIPOL's successful experience with a HWP 
will inform the design of future, more sensitive, 
experiments designed to characterize B-mode po- 
larization of the CMB. 
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Appendix: Beam Mapping 

The raw data collected during the Jupiter ob- 
servations contained the following systematic er- 
rors: the HWP synchronous signal, low-frequency 
drifts, and transients. These spurious signals 
needed to be removed to make accurate beam 
maps. The magnitude of the Jupiter signal was 
similar to the magnitude of the HWP synchronous 
signal, so the algorithm described here was tai- 
lored for recovering signals that were much greater 
than the noise. 

To remove the low- frequency drifts, we first re- 
moved a preliminary estimate of the HWP syn- 
chronous signal. This fi was modeled in this algo- 
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rithm as, 

8 

n{t) = + ^»^) ( "'^(^) + ®" ) • (1^) 
?i=i 

The fit parameters, Xn, Yn, and 0„ were es- 
timated using the foUowing two step procedure. 
First, to determine the phase of each Fourier com- 
ponent of the HWP synchronous signal, 0„, the 
following cosine-wave model was fit to the raw 
data binned in the HWP angle domain: 

8 

H(/3) = ^Z„cos(n/3 + e„). (18) 

ri=l 

Here, Zn and 0„ are the fit parameters, and 0„ 
in EquationfTTlis equal to 8„ in EguationfTSl The 
bin size was set to 1 deg. Given the HWP rotation 
frequency, this bin size prevented neighboring time 
samples from falling into the same bin. During 
the binning procedure, low-frequency noise was 
rejected in each bin by high-pass filtering the raw 
data in the frequency domain before binning; tran- 
sients, and planet signals were rejected by setting 
the bin value equal to the mode, which was esti- 
mated by iteratively histogramming the bin data. 
In the second step of estimating the f5, the linearly 
varying amplitude of each Fourier component was 
found by fitting a line to the demodulated data in 
the time domain. For demodulation, the reference 
signal, 

7^(^) = cos ( n/3(i) + e„ ) , (19) 

was phase locked by construction because the best- 
fit phase, which was output from the previous step, 
was used for each Fourier component. An esti- 
mate of the HWP synchronous signal was then 
constructed using 0„ from step one and X„ and 
Yn from step two. This estimate was then sub- 
tracted from the raw data leaving only Jupiter 
signal, transients, low-frequency drifts, and noise. 

Low- frequency drifts, which biased the beam 
maps if not subtracted, were removed from this 
data by iteratively fitting and subtracting second- 
degree polynomials from 24 sec long segments of 
data. This segment length was selected because it 
was much longer than a typical crossing time of 
Jupiter through the beam. After each iteration, 
the polynomial fit was subtracted and data points 
greater than 3 standard deviations away from zero 
were ignored in the subsequent fitting iterations. 



Given the 3 standard deviation rejection criteria, 
the polynomial estimates converged after three it- 
erations. This masking procedure prevented tran- 
sients and the planet signal from biasing the esti- 
mate of the low-frequency drifts. 

The final drift estimate was subtracted from the 
raw data. The effects of the electronic filters and 
the bolometer time constants were deconvolved us- 
ing the procedure given in Section 15.21 A second 
iteration of HWP synchronous signal estimation 
was required because the fi was phase shifted by 
the filter deconvolution. Transients larger than 
the empirically determined maximum Jupiter sig- 
nal were flagged. The beams were then mapped 
using the telescope pointing, the transient flags, 
and the data with fi and low-frequency drifts sub- 
tracted. The map pixel size was set to 0.7 arcmin 
to allow accurate estimation of the Bi, which was 
used during CMB power spectrum estimation. 

Biases in the beam map that were introduced 
by the low-frequency drifts were significant. To 
allowed for better estimation and subtraction of 
these drifts, the entire process outlined above was 
repeated with a Jupiter signal template removed 
during the low-frequency drift estimation. This 
signal template was computed by scanning a beam 
model with the telescope pointing. Here the beam 
model was the two-dimensional elliptical Gaussian 
that best fit the pixelized maps that were output 
from the first iteration of the process. 
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